Balance network of asymmetric simple exclusion process 
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We investigate a balance network of the asymmetric simple exclusion process (ASEP). Subsystems 
consisting of ASEPs are connected by bidirectional links with each other, which results in balance 
between every pair of subsystems. The network includes some specific important cases discussed 
in earlier works such as the ASEP with the Langmuir kinetics, multiple lanes and finite reservoirs. 
Probability distributions of particles in the steady state are exactly given in factorized forms ac- 
cording to their balance properties. Although the system has nonequilibrium parts, the expressions 
are well described in a framework of statistical mechanics based on equilibrium states. Moreover, 
the overall argument does not depend on the network structures, and the knowledge obtained in 
this work is applicable to a broad range of problems. 
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I. INTRODUCTION 

The asymmetric simple exclusion process (ASEP) is 
one of the most paradigmatic models to understand phe- 
nomena in nonequilibrium physics The model, con- 
sisting of a one-dimensional lattice and particles with 
hard-core exclusion interaction, describes fundamental 
transport phenomena and is applied to a broad ranee 
of problems: traffic flow 0], biological transport 0, [B~ 
0], and etc. As natural extensions of the ASEP, the 
effects of particle attachments and detachments in the 
bulk HHU, and multiple lanes [M3 have been inves- 
tigated, and some significant results have been presented. 
These systems allow additional motion of particles in the 
ASEPs and can be interpreted as networks of the ASEPs 
and reservoirs, where each site in the lattice is connected 
with the particle reservoir or a site in the different ASEP. 
On the other hand, the ASEP on networks has been fo- 
cused on recently |16l - [l9| . The results have concluded 
that the dynamics of the system depends on structure 
of the networks. In this paper, we focus on an exactly 
solvable network consisting of the periodic ASEPs. The 
steady state of the system is described by general expres- 
sions, which are found to be independent of the topology 
of the network. The key to construct the expressions 
is the detailed balance satisfied among the subsystems, 
namely, the ASEPs are in balance with each other in the 
network. This kind of structure has been reported in the 
previous studies @, Q , and we have successfully general- 
ized the system in this work. We provide a certain class 
of solvable TASEP systems, paving the way for the use 
of this structure. 

The rest of this paper is organized as follows. Section 
[H] gives the definition of the network and the relationship 
to its related models. In Sec. IIII1 we give the exact sta- 
tionary distribution of the system. Using the expressions, 



we derive some physical quantities in Sec. IIVI Finally, we 
summarize the discussion in Sec. IVl 
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FIG. 1. (color online). Examples of the balance network. The 
balance network can be regarded as generalization of (b) the 
ASEP with the Langmuir kinetics, (c) the multilane ASEP, 
and (d) the simple exclusion network. 



II. MODEL 

We consider a network of exclusion process consisting 
of particles, sites, links and a single reservoir. Each site 
can contain at most one particle, and each particle jumps 
to a site in the same subsystem or to a site in another 
subsystem through a link. Here the ASEP on a ring is 
mainly focused on as the subsystem. The periodic ASEP 
has an 'equiprobable' property that all the configurations 
of particles, {t?}j, appear equally likely in the steady 
state, provided the number of particles is fixed. Here, 
t\ is the occupation number of site i (i = 1, • • • , Lj) in 
the subsystem j (j — 1, ■■■,!(), and {r-}j is a set of 
the occupation numbers that describes each configura- 
tion in subsystem j. In principle, other exclusion pro- 
cesses (or even processes of bosons discussed later) can 
also be candidates for this subsystem if only they satisfy 
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the equiprobable property. In this work, we also con- 
sider a single site without dynamics in itself as one of the 
equiprobable subsystems. The components in the system 
are summarized as follows: 

(i) Equiprobable subsystem: A set of sites that has 
equiprobable dynamics such as the ASEP with periodic 
boundary conditions. Each site in an equiprobable sys- 
tem j has a common leaving rate of particles, Xj'i 

(ii) Link: Bidirectional links connect pairs of sites in 
different subsystems (unidirectional links are forbidden). 
Each site can have an arbitrary number of links; 

(iii) Reservoir: A reservoir can accept and provide an 
arbitrary number of particles through links, and its pro- 
vision rate is XR- Only a single reservoir is allowed in the 
system. 

These components are set in the system as a network 
[see Fig. [1] (a)] . Note that the network must be the con- 
nected network: the network cannot have isolated parts, 
and every pair of subsystems must be interconnected by 
links and/or other subsystems. In the network particles 
jump to the neighboring sites, following the hard-core ex- 
clusion principle. As shown in Fig. [5] (a), a particle at 
a randomly chosen site in the subsystem j jumps to the 
next site to its right with a rate pj, and to the linked 
sites in other subsystems with a rate \ji if the target 
sites are empty. Moreover, through the links, the reser- 
voir can accept and provide particles with rates Xj an d 
Xr, respectively [see Fig. H](b)]. 

By these formulations of the system, we can see that 
the network includes some important cases; the ASEP 
with the Langmuir kinetics [3] [Fig. EDO)], the multilane 
ASEP d, [To[ [Fig. Q] (c)] , and a simple exclusion network 
[Fig. [2(d)]. Furthermore, the balance network generally 
represents the multiple competing ASEPs. In the context 
of biology, the competition of the ASEPs is discussed as 
the problem of multiple mRNAs [20], or it may explain 
the dynamics of motor proteins on a spindle consisting 
of microtubules in cell division. 
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FIG. 2. (color online). Transition rules of the system. Parti- 
cles hop to the next site with a rate pj (the ASEP) and leave 
the subsystem j with a rate Xi through the links, obeying the 
exclusion principle. The reservoir contains an infinite number 
of particles and provides a particle to an empty linked site 
with a rate xr- 



The probability of finding the system in a configuration 
({Tih, ■ • • , {t?}k) is given by 



III. EXACT ANALYSES 

We analyze the balance network in the steady state, 
focusing on the probability distribution for each config- 
uration of particles. First let us review the expressions 
for the ASEP on a ring. A possible configuration {Ti} is 
realized with the probability 
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P({n})=N- 1 f({r i }) ) 



(1) 



where /({r^}) is the probability weight of each configura- 
tion, and N -1 is the normalization factor. Since all the 
possible configurations in this system are equally likely in 
this system, /({Tj}) = 1. This property does not depend 
on the system length, Lj, and the density of particles. 

Then, we present the probability distributions for the 
balance network using the weight of each configura- 
tion in subsystem j with rij particles, /n, ({ T f}j) = 1- 
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Here S _1 is the normalization factor. As shown in the 
next part, this S corresponds to the grand partition func- 
tion in statistical mechanics. Note that, even if the sys- 
tem lacks the reservoir, this expression can be used with 
a slight modification. In this case, since the absence of 
the reservoir leads to a constraint on the particle number, 
the sum of the weights is taken over all the configurations 
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with a given number of particles, n, while 2 is obtained 
by considering all the possible configurations for any par- 
ticle numbers. Furthermore, the provision rate \R does 
not influence the equations because each Xr in Eq. © 
is cancelled out by the normalization factor. Although 
the conservative systems are also interesting when we 
consider actual biological processes with finite resource 
[t| |20| - |22| . constraints on particle numbers often cause 
difficulty of analysis. 

Let us confirm that these expressions correctly describe 
the system in the steady state by considering the master 
equation: 



= T P(C) = J2 {P(C')W(C -»• C) - P(C)W(C -¥ C')}, 

(4) 

where C and W(C — s> C) indicates the configuration of 
particles and the transition probability from configura- 
tion C to C, respectively. Here we separate the transi- 
tions into three parts, i.e., internal transitions in each 
subsystem, intersubsystem transitions, and transitions 
between the reservoir and subsystems. Since each in- 
ternal transition does not change the particle numbers 
in the subsystems, it is obvious that Eq. ^ satisfies the 
master equation for these transitions (for each subsystem, 
the equiprobable expression for a fixed particle number 
satisfies the master equation for the equiprobable sub- 
system, and thus the internal transition terms vanish). 
Generally, these terms vanish through taking the sum of 
all the transitions, and the detailed balance conditions 
are not satisfied: this is the generalization of the detailed 
balance to the nonequilibrium steady state [Hj]. On the 
other hand, the other two types of transitions satisfy the 
detailed balance conditions: 



= P(C')W(C -»■ C) - P(C)W(C -> C). 



(5) 



Let us take a transition between subsystem j\ and j'2 
(ji < 32) as an example. In the transition, a par- 
ticle jumps from site i\ in subsystem j\ to site ii 
in the subsystem ji through a link, which results in 
a change of the particle numbers in the subsystems, 
{•••_> n h ,--- ,n h ,---} ->■ {■■■,n h -l,---,n h + l,---}. 
Taking its reverse transition into account, Eq. §5§ holds: 
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Each bidirectional link ensures the existence of the re- 
verse transition for a given intersubsystem transition. In 
the same manner, we can prove these detailed balance 
conditions for the transitions between the equiprobable 
subsystem and the reservoir. Moreover, these cancella- 
tion mechanisms are independent of the network struc- 
ture and capacity of sites; the expressions are valid for 



finite pools of particles 24] , and there is scope for exten- 
sion of the subsystems to multiple occupation processes. 

To summarize, all the terms in Eq. (j4j vanish accord- 
ing to the properties of transitions, i.e., the nonequilib- 
rium in the internal transitions and the balance in the 
external transitions. 

From Eqs. @ and we can derive the current of 
particles in subsystem j defined as Jj — (rj (1 — t/ + i)) 
by considering all the configurations with t\ — 1 and 
t' =0: 
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In the first expression we swapped the summation and 
the multiplication. The current is the quantity of great 
importance to describe transportation phenomena and 
the characteristic quantity in the nonequilibrium sys- 
tems. It is noteworthy that the current is determined 
only by the parameters of the subsystem and the reser- 
voir. Moreover, we can prove that the correlation be- 
tween the occupation numbers of successive two sites can 
be ignored even for finite size of the systems [25| . 



IV. CORRESPONDENCE TO STATISTICAL 
MECHANICS 

The structure of the Eq. (|2|) is well explained in a 
framework of statistical mechanics. Let us derive the 
expected value of the occupation numbers, putting Xj — 
e l3ej and \R = to emphasize the correspondence. The 
grand partition function is calculated as 

^^^■•En^ti 3 ') (9) 

n 1= n K =0j=l \ 3 / 
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where Ej is defined as 3j = (1 + e~P^ j ~^ ) Lj . Then, the 
expected value of the occupation number, (rij) is given 
by 

^■)=^E---E^n^'-^ 
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Thus, the density of particles in each subsystem is de- 
rived. Note that these calculations can be performed 
without the interpretation using the energy, the inverse 
temperature, and the chemical potential; however, the 
expressions are highly suggestive. If one regards each site 
as a distinctive energy state with energy, ej, of fermions, 
Eq. (|TB)) coincides with the Fermi distribution [let Lj = 1 
for the simplicity of the argument (Lj corresponds to the 
degeneracy)]. On the other hand, the system can also 
be interpreted as a problem of chemical adsorption with 
chemical potential, /i, and stabilization energy, —ej. In 



this case, Eq. (|16[) corresponds to the Langmuir isotherm 
of the system with independent Lj sites in contact with 
the reservoir. Since each pair of connected subsystems 
are in balance, the network is equivalent to a set of sep- 
arated subsystems in balance with the reservoir. Thus, 
the steady state is determined only by the parameters of 
each subsystem and the common reservoir. This is the 
reason why the overall argument can be well understood 
in the framework of statistical mechanics. However, it 
is still noteworthy that the statistical mechanics expres- 
sions can be naturally extended to the system consisting 
of some nonequihbrium parts. 

Additionally, we evaluate the variance of the particle 
number in the subsystem j, Var[ physical quan- 

tity characterizing the equilibrium steady state. 
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Interestingly, the expression is associated with the cur- 
rent Eq. © as 



Var[nj] = LjJ . 



(22) 



Thus, an intriguing relation is derived, where the rep- 
resentative quantities of equilibrium and nonequihbrium 
physics are linked together. [26| 



V. CONCLUSIONS 

We have presented the balance network consisting of 
nonequihbrium subsystems, bidirectional links, and a sin- 
gle reservoir. The network includes a wide variety of 
models relevant to previous works and is very useful in 
the meaning of application. On the other hand, the net- 
work has a prominent structure of balance connections, 
which allows us analytical solutions. From the probabil- 
ity distribution of particles we can calculate some phys- 
ical quantities, and the overall argument can be well 
understood in the framework of established statistical 
physics. 

In the balance network, only the bidirectional links are 
allowed because unidirectional links will cause the 'flow' 
of particles between subsystems and violate the balance 
relations. Besides, if we allow a single site which can 
contain more than one particle, the site is equivalent to 
a finite pool of particles or a finite reservoir. The bal- 
ance network can contain an arbitrary number of finite 
reservoirs; on the other hand, construction of exact prob- 
ability distribution for the system with multiple infinite 
reservoirs is not straightforward. In future works, further 
analyses on the extension of the balance network and its 
relations with nonequihbrium physics are needed. 
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